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Abstract: On Barabasi-Albert networks with z neighbours selected by each 
added site, the Ising model was seen to show a spontaneous magnetisation. This 
spontaneous magnetisation was found below a critical temperature which increases 
logarithmically with system size. On these networks the majority-vote model with 
noise is now studied through Monte Carlo simulations. However, in this model, the 
order-disorder phase transition of the order parameter is well defined in this system 
and this wasn't found to increase logarithmically with system size. We calculate 
the value of the critical noise parameter q c for several values of connectivity z of 
the undirected Barabasi-Albert network. The critical exponentes (3/v, j/v and 
1/v were calculated for several values of z. 
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Introduction 

It has been argued that nonequilibrium stochastic spin systems on reg- 
ular square lattice with up-down symmetry fall in the universality class of 
the equilibrium Ising model pQ. This conjecture was found in several models 
that do not obey detailed balance [21 El El- Campos et al. |5J investigated the 
majority- vote model on small-world network by rewiring the two-dimensional 
square lattice. These small-world networks, aside from presenting quenched 
disorder, also possess long-range interactions. They found that the critical 
exponents j/v and /3/v are different from the Ising model and depend on 
the rewiring probability. However, it was not evident whether the exponent 
change was due to the disordered nature of the network or due to the pres- 
ence of long-range interactions. Lima et al. jH] studied the majority-vote 
model on Voronoi-Delaunay random lattices with periodic boundary condi- 
tions. These lattices posses natural quenched disorder in their conecctions. 
They showed that presence of quenched connectivity disorder is enough to 
alter the exponents (3/v and 7/1/ the pure model and therefore that is a 
relevant term to such non-equilibrium phase-transition. Sumour and Shabat 
0|H] investigated Ising models on directed Barabasi-Albert networks with 



1 



the usual Glauber dynamics. No spontaneous magnetisation was found, in 
contrast to the case of undirected Barabasi-Albert networks PUJ El Ej where 
a spontaneous magnetisation was found lower a critical temperature which 
increases logarithmically with system size. Lima and Stauffer ^3] simulated 
directed square, cubic and hypercubic lattices in two to five dimensions with 
heat bath dynamics in order to separate the network effects form the effects 
of directedness. They also compared different spin flip algorithms, including 
cluster flips > for Ising-Barabasi- Albert networks. They found a freezing- 
in of the magnetisation similar to [JJ |H] , following an Arrhenius law at least 
in low dimensions. This lack of a spontaneous magnetisation (in the usual 
sense) is consistent with the fact that if on a directed lattice a spin Sj in- 
fluences spin Si, then spin 5, in turn does not influence Sj, and there may 
be no well-defined total energy. Thus, they show that for the same scale-free 
networks, different algorithms give different results. More recently, Lima [Tf] 
investigated the majority-vote model on directed Barabasi-Albert network 
and calculated the j3/v, 7/V, and 1/v exponents and these was different 
from the Ising model and depend on the values of connectivity z of the di- 
rected Barabasi-Albert network. Now, we calculate the same (3/v, 7/V, and 
1/v exponents for majority- vote model on undirected Barabasi-Albert net- 
work and these are different from the Ising model and depend on the values 
of connectivity z of the directed Barabasi-Albert network. 



Model and Simulaton 

We consider the majority-vote model, on directed Barabasi-Albert Net- 
works, defined [151117101111 by a set of "voters" or spins variables a taking 
the values +1 or —1, situated on every site of an undirected Barabasi-Albert 
Networks with N sites, and evolving in time by single spin-flip like dynamics 
with a probability Wi given by 
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where S(x) is the sign ±1 of x if x 7^ 0, S(x) = if x = 0, and the sum runs 
over all nearest neighbours of <7j. In this network, each new site added to the 
network selecting z, already existing sites as neighbours influencing it; the 
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newly added spin does not influence these neighbours. The control parameter 
q plays the role of the temperature in equilibrium systems and measures the 
probability of aligning antiparallel to the majority of neighbours. 

To study the critical behavior of the model we define the variable m = 
Y,iLi a i/N. In particular , we were interested in the magnetisation, suscep- 
tibility and the reduced fourth-order cumulant: 

M(q) = [< \m\ >U (2) 

X (q) = N[<m 2 >-< \m\ >\ V) (3) 



U(q) 



< m 4 > 



1 

av 



(4) 



3 < \m\ > 2 

where < ... > stands for a thermodynamics average and [...] a „ square brackets 
for a averages over the 20 realizations. 

These quantities are functions of the noise parameter q and obey the 
finite-size scaling relations 

M = N-^f m (x)[l + ...], (5) 
X = N^f x (x)[l + ...], (6) 

^ = N 1 '»f u (x)[l + ...], (7) 

where z/, (3, and 7 are the usual critical exponents, fi(x) are the finite size 
scaling functions with 

x=(q- q^N 1 ^ (8) 

being the scaling variable, and the brackets [1 + ...] indicate corretions-to- 
scaling terms. Therefore, from the size dependence of M and x we obtained 
the exponents (3/v and 7/z/, respectively. The maximum value of susceptibil- 
ity also scales as N 1 ^ . Moreover, the value of q for which x has a maximum, 
qXmax — q c (N), is expected to scale with the system size as 

q c (N) =q c + bN- 1 '", (9) 

were the constant b is close to unity. Therefore, the relations (7) and (9) are 
used to determine the exponente 1/z/. We have checked also if the calculated 
exponents satisfy the hyperscaling hypothesis 

2(3/u + 1 /u = D eff (10) 
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in order to get the effective dimensionality, -D e //> for various values of z. 

We have performed Monte Carlo simulation on undirected Barabasi- Albert 
networks with various values of connectivity z. For a given z, we used sys- 
tems of size N = 1000, 2000, 4000, 8000, and 16000. We waited 10000 Monte 
Carlo steps (MCS) to make the system reach the steady state, and the time 
averages were estimated from the next 10000 MCS. In our simulations, one 
MCS is accomplished after all the N spins are updated. For all sets of pa- 
rameters, we have generated 20 distinct networks, and have simulated 20 
independent runs for each distinct network. 

Results and Discussion 

In Fig. 1 we show the dependence of the magnetisation M and the sus- 
ceptiblity x on the noise parameter, obtained from simulations on undirected 
Barabasi- Albert network with 16000 sites and several values of connectivity 
z. In the part (a) each curve for M, for a given value of N and z, suggests 
that there is a phase transition from an ordered state to a disordered state. 
The phase transition occurs at a value of the critical noise parameter q c , 
which is an increasing function the connectivity z of the directed Barabasi- 
Albert network. In the part (b) we show the corresponding behavior of the 
susceptibility x, the value of q where \ has a maximum is here identified as 
q c . In Fig. 2 we plot Binder's fourth-order cumulant for different values of 
N and two different values of z. The critical noise parameter q c , for a given 
value of z, is estimated as the point where the curves for different system 
sizes N intercept each other. In Fig 3 the phase diagram is shown as a func- 
tion of the critical noise parameter q c on connectivity z obtained from the 
data of Fig. 2. 

The phase diagram of the majority-vote model on undirected Barabasi- 
Albert network shows that for a given network (fixed z ) the system becomes 
ordered for q < q c , whereas it has zero magnetisation for q > q c . We notice 
that the increase of q c in function the z is slower of the one than in ^H]. In the 
Fig. 4 and 5 we plot the dependence of the magnetisation and susceptibility, 
respectively at q = q c with the system size. The slopes of curves correspond to 
the exponent ratio f3/v and ■y/u of according to Eq. (5) and (6), respectively. 
The results show that the exponent ratio (3/v increase and 7/V decrease at 
q c when z increase, see Table I. 

In Fig. 6 we display the scalings for susceptibility at q = q c (N) (square), 
x{<lc{N)), and for its maximum amplitude, XN aX i an d the scalings for suscep- 
tibility at the q = q c obtained from Binder's cumulant, x(<?c) (circle), versus 



4 



N for connectivity z = 8. The exponents ratio 7/z/ are obtained from the 
slopes of the straight lines. For almost all the values of z, the exponents 7/z/ 
of the two estimates agree (along with errors), see Table I. An increased z 
means a tendency to decrease the exponent ratio 7/z/, see Table I, they agree 
with the results of Luiz et al [THj . but disagree with the results of Lima 14: 
for directed Barabasi- Albert network, where the values of the exponents ratio 
7/z/ are all different and with a slight tendency to increase at q c and decrease 
at q c {N). Therefore we can use the Eq. (9) , for fixed z, obtain the critical 
exponent 1/z/, see Fig. 7. In the Fig. 8 we show the critical behavior of /3/u, 
7/z/ and 1/z/ the exponentes in function of connectivity z. 

To obtain the critical exponent 1/z/, we calculated numerically U (q) = 
dU(q)/dq at the critical point for each values of iV at connectivity fixed z. 
The results are well in agreement with the scaling relation (7) . Then , we 
also can calculate the exponents 1/z/, through this relation. Therefore we do 
not have to get the values of the exponents 1/z/ for each connectivity z 

The Table I summarizes the values of q c , the exponents (3/v, 7/z/, and the 
effective dimensionality of systems. For all values of z, D e ff — 1 5 which has 
been obtained from the Eq. (9), therefore when z increases, (3/v increases 
and 7/z/ decreases at q c , thus providing the value of D e ff = 1 (along with 
errors). Therefore, the undirected Barabasi-Albert network has the same 
effective dimensionality that Erdos-Renyi's random graphs [1§| and directed 
Barabasi-Albert network ^3]. J. M. Oliveira ^H] showed which majority- 
vote model defined on regular lattice has critical exponents that fall into 
the same class of universality as the corresponding equilibrium Ising model. 
Campos et al [3] investigated the critical behavior of the majority- vote on 
small-world networks by rewiring the two-dimensional square lattice, Luiz et 
al [TH| studied this model on Erdos-Renyi's random graphs, and Lima et al 
[U] also studied this model on Voronoy-Delaunay lattice and Lima directed 
Barabasi-Albert network The results obtained these authors show that 
the critical exponents of majority-vote model belong to different universality 
classes. 

Finally, we remark that our MC results obtained on undirected Barabasi- 
Albert network for majority-vote model show that critical exponents are 
different from the results of JHI for regular lattice, of Luiz et al [TS] for 
Erdos-Renyi's random graphs and Lima [T4*] . 

Conclusion 

In conclusion, we have presented a very simple nonequilibrium model on 
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3 


0.259(2) 


0.133(21) 


0.713(18) 


0.655(31) 


0.83(7) 


0.979(27) 


4 


0.306(3) 


0.231(22) 


0.537(8) 


0.519(17) 


0.43(2) 


0.999(23) 


6 


0.355(2) 


0.283(8) 


0.445(15) 


0.423(3) 


0.35(5) 


1.011(17) 


8 


0.380(6) 


0.323(2) 


0.358(7) 


0.405(6) 


0.39(3) 


1.004(7) 


10 


0.396(3) 


0.338(2) 


0.324(2) 


0.380(3) 


0.324(5) 


1.000(2) 


20 


0.428(2) 


0.344(2) 


0.305(2) 


0.350(2) 


0.307(5) 


0.993(2) 


50 


0.456(3) 


0.366(2) 


0.255(2) 


0.341(3) 


0.30(1) 


0.987(2) 


100 


0.471(3) 


0.373(2) 


0.218(5) 


0.330(3) 


0.308(4) 


0.964(5) 



Table 1: The critical noise q c , the critical exponents, and the effective dimen- 
sionality D e ff , for undirected Barabasi- Albert network with connectivity z. 
Error bars are statistical only. 

undirected Barabasi- Albert network [7J Ej. Different of of Ising model, in 
these networks, the majority-vote model presents a second-order phase tran- 
sition which occurs with connectivity z > 1. The exponents obtained are 
different from the other models. Nevertheless, our Monte Carlo simulations 
have demonstrated that the effective dimensionality D e ff equals unity, for 
all values of z, which agrees with the results de Luiz et al JH]- However, 
when z grows, the exponents in the critical point q c , (3/v obtained through 
Binder's cumulant grows and exponents 7/1/ decrease, satisfying the hyper- 
scaling relation with D e ff = 1- 
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n=16000 sites 




Figure 1: Magnetisation and susceptibility as a function of the noise param- 
eter q, for N = 16000 sites. From left to right we have z — 2, 3, 4, 6, 8, 10, 
20, and 50 . 
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Figure 2: Binder's fourth-order cumulant as a function of q. We have z = 2 
and z = 50 for iV = 1000, 2000, 4000, 8000 and 16000 sites. 
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Figure 3: The phase diagram, showing the dependence of critical noise 
rameter q c on connectivity z. 
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Figure 4: In M(q c ) versus In N. From top to bottom, z — 2, 3, 4, 6, 8, 10, 
20, and 50. 
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Figure 7: Plot of In \q c (N) — q c \ versus In N. From bottom to top z — 2, 3, 
4, 6, and 8. 



14 



1.0 
0.9 
0.8 
0.7 
> 0.6 
-> 0.5 

0.3 
0.2 
0.1 
0.0 





o-- 


- - G 


, P/v 


QX 

- \ '\ 
{ i 


B - - 


- - Q 


, y/v 


- Xn 

. 

II 
H 
1 


x- - - 


- - X 


, 1/v 


1 

II 
11 

ED 

i \ 








■ *\ 














-o-- 




P 

©' 

/ 






-E 


/ 

/ 

© 

/ 

/ 








/ 

o 

1 , 1 









10 



20 



30 
z 



40 



50 



60 



Figure 8: Critical behavior the f3/v, 7///, and 1/v exponents in function of 
connectivity z. 



15 



